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The “d-wave” symmetry of the supereondueting order 
in the euprate high temperature supereonduetors is a well 
established faet [1, 2], and one whieh identifies them 
as ^'uneonventional” However, in maeroseopie eontexts 
- ineluding many potential applieations (i.e. supercon- 
dueting ^^wires^’) - the material is a eomposite of ran¬ 
domly oriented supereondueting grains in a metallie ma¬ 
trix, in whieh Josephson eoupling between grains mediates 
the onset of long-range phase eoherenee. (See, e.g., [3- 
5].) Here, we analyze the physies at length seales large 
eompared to the size of sueh grains, and in partieular 
the maeroseopie eharaeter of the long-range order that 
emerges. While XY-glass order and maeroseopie d-wave 
supereonduetivity may be possible, we show that under 
many eireumstanees - espeeially when the d-wave super¬ 
eondueting grains are embedded in a metallie matrix - 
the most likely order has global s-wave symmetry. 
Classification of phases: 

The anomalous average of the spin-singlet electron 
pair-creation operator which characterizes the supercon¬ 
ducting state is [? ] 

r)) = + ip^{r')ipi{r)])/^2, ( 1 ) 

where (r) creates an electron with spin polarization a 
at position r, and ( ) represents the equilibrium average 
over all thermal and quantum fluctuations. 

Above the critical temperature, Tc, in the normal 
metal phase, ( 0 (r, f^)) = 0 , while ( 0 (r, r^)) 7 ^ 0 for all 
T <Tc. In pure crystals this quantity is only a function 
of {r — r'), so it is convenient to introduce its Fourier 
transform {(j){p)). Possible superconducting phases were 
classified in Refs. [ 6 , 7]. In particular, in s-wave su¬ 
perconductors, ( 0 (p)) is invariant under all symmetry 
transformations of the crystal and ( 0 (r, r)) 7 ^ 0 , while 
for other forms of singlet order, ( 0 (p)) changes sign un¬ 
der certain symmetry transformations, and consequently 
(0(r, r)) = 0. For example, in a d-wave superconductor, 
( 0 (p)) (as well as {(j){r,r'))) changes sign under rotation 
by 7 r/ 2 , as is shown schematically in Fig. 2. 

In a disordered system, {(j){r,r')) is a sample specific 
random quantity which does not posses any spacial sym¬ 
metry. As a result, in any singlet superconducting phase 
in a disordered system, since no symmetry prevents it, 
generically ( 0 (r, r)) 7 ^ 0 . [? ] 

In this article we assume that the disorder ensemble 
preserves a group of translational and rotational symme¬ 
tries statistically {i.e. on average). We will show below 
that after averaging over the realizations of the disor¬ 


der, several superconducting phases emerge that can be 
precisely characterized by different order parameters: 

The internal symmetry of the superconducting pairs 
(for example, s-wave or d-wave) refers to the transforma¬ 
tion properties of the quantity 

^{r — H) = {(j){r,r')). ( 2 ) 

under rotation, where the overline indicates a quan¬ 
tity that has been averaged over configurations of the 
quenched variables, e.g. the size, shape, orientation, and 
location of the superconducting grains as well as over 
realizations of the disordered potential in the metal. 

Since the internal structure of ^(r) is not directly mea¬ 
surable, the symmetry of the superconducting state is 
best defined in terms of “phase sensitive measurements.” 
Deflnitionally, this refers to any measurement of the rela¬ 
tive phase of the order parameter at two macroscopically 
separated locations on the surface of the system. For in¬ 
stance, this can be measured in a SQUID consisting of 
Josephson junctions at two surface positions connected 
by an external (macroscopic) conventional superconduct¬ 
ing wire loop. In an s-wave state, the phase difference 
around the loop is zero independent of the relative ori¬ 
entation of the two surfaces; consequently in equilibrium 
there is no magnetic flux through the SQUID. By con¬ 
trast, in a d-wave state, a phase difference of tt is induced 
in a “corner SQUID” in which the loop connects portions 
of the surface approximately at right angles to each other, 
so there is a half quantum of magnetic flux through the 
SQUID in equilibrium [I, 2 ]. At the same time, if the 
SQUID loop connects nearly parallel portions of surface 
of a d-wave state (either on the same or opposite sides 
of the system), no equilibrium flux is induced. Similar 
analysis can be used to give a phase sensitive definition 
of other possible pairing symmetries. 

The appropriate order parameter characterization of a 
superconducting glass state is a bit more subtle: while 
{(j){r,r')) 7 ^ 0 , its local phase varies randomly as a func¬ 
tion of position, and correspondingly its configuration (or 
spatial) average vanishes. However, this state is sharply 
distinguished from the normal state by the existence of 
a non-zero Edwards-Anderson type order parameter [ 8 ], 

M(f-f') = |(</>(f,F))|2, (3) 

i.e. the glass state has {(l){f,r')) = 0 but M = M(0) > 0. 

Another feature of a glassy state, reflecting the exis¬ 
tence of random variations of the phase of ( 0 (r, U)), is 
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that there are local equilibrium currents, {J{r)) 7 ^ 0 , and 
associated spontaneous breaking of time-reversal symme¬ 
try. It is easy to see that the configuration (or spatial) 

average current must vanish, {J{r)) = 0. Instead, an 
appropriate tensor order parameter is 

Tah = (Ja(0))(J6(0)). (4) 

A complete classification of the various broken symme¬ 
try phases in disordered superconductors is not currently 
available. However, in terms of the various quantities in¬ 
troduced above, we can define the phases which we will 
encounter in our discussion: 

• 1) In the “normal” state at elevated temperatures, 
no symmetries are broken, and hence all the order 
parameters vanish. 

• 2) In a “s-wave superconducting” state, 4>(r) is non¬ 
zero, while Ta^h = 0 - z.e. there are no equilibrium 
currents in the bulk, and 4>(r) is invariant under 
all symmetry transformations, i.e. it is rotation- 
ally invariant. In phase sensitive measurements, no 
equilibrium flux is induced in a SQUID, regardless 
of its geometry. 

• 3) In a “d-wave superconducting” state, Ta^b = 0, 
4>(r) is non-zero for non-zero |r| and changes sign 
under rotation by 7r/2 about some axis, but is in¬ 
variant under rotation by tt. Consequently, 4>(0) = 
0 and a half flux quantum is induced in equilibrium 
in a suitable corner SQUID. 

• 4) In a “superconducting XY glass” state, 4>(r) van¬ 
ishes but both M and Ta^b are non-zero. It is widely 
accepted, but still not completely settled that such 
a state exists in three dimensions below a non¬ 
zero glass transition temperature. There is also the 
possibility of a partially ordered glass phase with 
M(r) = 0 but Ta^b 7 ^ 0 , but whether this arises in 
generic models is still being debated. [9, 10] 

• 5) It is possible to have a phase of coexisting uni¬ 
form and glassy order - in such a phase both 4> 
and Tab are non-zero. Depending on the behav¬ 
ior of 4>(r) under symmetry transformations (“ro¬ 
tations”), such a phase can still be classified as s- 
wave or d-wave etc. 

• 2b) The most unexpected new state we identify 
here is what we will henceforth refer to as a ’‘’‘glob¬ 
ally s-wave superconducting’’^ state. Definitionally, 
this is simply an unfamiliar limit of a s-wave state 
- one in which locally (i.e. in each superconducting 
grain) the order parameter has d-wave symmetry, 
but globally it is s-wave. Spectroscopically, such 
a state can reflect its microscopic origins as a d- 
wave superconductor, but from the viewpoint of 


macroscopic phase-sensitive measurements, it has 
s-wave symmetry. Moreover, time-reversal symme¬ 
try is unbroken. Tab = 0. 

The effective Hamiltonian: Below the bulk transi¬ 
tion temperature, in a system composed of grains of size 
large compared to the coherence length, there is a well 
developed magnitude of the order parameter on each 
grain. We will consider the case of most relevance to 
the cuprates, in which the order parameter on each grain 
transforms according to a one-dimensional non-trivial (d- 
wave) representation of the point group. The only im¬ 
portant low energy degree of freedom is the overall phase 
of the superconducting order parameter on each grain, 
which we will designate Oj. 

There is a degree of arbitrariness in the definition of 
Oj; we choose a convention such that when Oj = 0, the 
order parameter on the grain is real and is positive in 
some particular crystalographic direction. Moreover we 
assume that the grains are sufficiently large that we can 
neglect quantum fluctuations of the order parameter. In 
this case, the macroscopic properties of the system can 
be captured by the phenomenological model 

H = -Y,Jijcos{8i-0j) , (5) 

ij 

where the Josephson coupling Jij between grains i and j 
is real. 

On general phenomenological grounds, we can separate 
the contributions to Jij into three pieces, 

(6) 

Here, > 0 for all a, while rjj = ±1 and rjij = ±1 are 
random variables with vanishing mean which determine 
the sign of the corresponding contributions to Jij. The 
magnitude of each contribution to Jij is characterized 

by its mean, Jif^ ^ and to be explicit, we will 

always consider this model in three spatial dimensions, 
although it need not be isotropic. 

For grains of conventional superconductors, under 
most circumstances, the Josephson coupling would be 
positive, which is to say that ~ ~ 0; in this 

limit, the model is equivalent to an XY ferromagnet 
which has a single phase transition at Tc ^ J^^\ where 
physically the “ferromagnetic phase” corresponds to a 
statistically uniform s-wave superconducting phase. (A 
mean-field estimate yields = J^^^/2.) 

In any case in which the crystalline axes of d-wave 
superconducting grains are embedded into a disordered 
metal with random orientations, the sign of Jij is ran¬ 
dom, which means that = 0. The two remaining 
terms in Eq 6 have very different character: The term 
proportional to J^j^ has its sign determined by a prod¬ 
uct of quantities that depend on the properties of each 
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grain separately (which we will see is roughly related to 
the shape of the grains). Conversely, the sign of the term 
proportional to Jh ^ is determined by a joint property of 
the pair of grains (which we will see is related to the rel¬ 
ative orientation of their crystalline axes). In the limit 
= 0 (and = 0), this problem is a version of 
the standard model of an XY spin-glass, [17] while for 
j(2) = 0 (and = 0), this problem is a version of the 
well known Mattis model. [11] 

Let us redefine the zero of phase on each grain sepa¬ 
rately to introduce “Mattis-transformed” phases 

0j=0j+Ti{l-r]j)l2. ( 7 ) 

In terms of these transformed variables, the form of H 
is unchanged, but the role of the different distributions 
is interchanged such that ^ ^ J^‘^\ and 

That the pure Mattis model is transformed 
in this way into a pure ferromagnetic XY model reflects 
the well known fact that this model introduces disorder 
without frustration. 

Conjectured phase diagrams: In Fig. la we show a 
conjectured phase diagram for the model in Eq. 5 under 
the conditions that = 0 as a function of the dimen¬ 
sionless temperature, T/ J, and the relative magnitude of 
the Mattis and spin-glass type couplings, 0 < /T < 1, 
where J = ^[JC)]2 + [J(2)]2, We have labeled the or¬ 
dered state at small / J “globally s-wave” as in the 

present context, as we will show below, this state - in 
which the Mattis transformed phases are uniformly (“fer- 
romagnetically”) ordered - corresponds to the state de¬ 
fined in 2b), above. The “XY glass phase” is defined in 
4), above. The intermediate state, in the present context, 
has coexsiting global s-wave and XY glass order, corre¬ 
sponding to 5), above. Formally, via the transformation 
in Eq. 7, the same phase diagram applies to the prob¬ 
lem (which has been studied in the spin-glass literature [? 
]) of an XY spin glass with an excess of ferromagnetic 
interactions (J^^^ = 0, while and ^ 0). 

The arguments leading to this phase diagram along the 
edges are as follows: 

i) For = 0 the Mattis transformed problem is 
equivalent to a XY ferromagnet with some randomness 
in the magnitude of the exchange couplings. Thus, we 
conclude that there is a single phase transition with 
Tc ^ to a phase with long-range superconducting 
order. 

ii) Decades of work has still not resulted in a well es¬ 
tablished understanding of even the most basic features 
of the XY spin-glass in d = 3. It is well accepted that 
there is a thermodynamic transition to a XY glass phase 
in large enough d. (See for example Ref. [9].) However, 
while it is widely believed that this conclusion applies 
in d = 3 (but not in d = 2), there remains some un¬ 
certainty concerning this conclusion. [10]. Numerical ex¬ 
periments certainly reveal that if there is a transition at 


Tgiass oc , the proportionality constant must be small 
(z.e. the spin-glass transition temperature is at least an 
order of magnitude smaller than for the corresponding 
model without frustration). Another intensely debated 
issue is whether there exists an intermediate partially or¬ 
dered phase with M = 0 but Tab ^ 0. We have drawn 
the phase boundaries in the limit ^ 0 in Figs, la 
and b under the assumption that there is a finite Tglass 
in d=3. If it turns out that T,lass = 0, the lines mark¬ 
ing the boundaries of the various glassy phases should be 
reinterpretted as crossover lines below which relaxation 
rates become extremely small. Conversely, if there are 
two transitions, then the lines should be interpretted as 
the mean of the two transitions. Even if there is a tran¬ 
sition, it is not clear whether the spin-glass phase would 
have a non-zero critical current. Moreover, in a glassy 
state, thermodynamic and transport quantities (includ¬ 
ing any “apparent critical current”) would be time de¬ 
pendent. This is very different from the ’’globally s-wave 
superconducting state” where the critical current is finite 
and time independent. Approaching the glass phase from 
above (T > Tgiass)^ the tendency to a state with non-zero 
Tab should give rise to a growing paramagnetic response 
- related to the so-called “Wohlleben effect.” [18, 19] 

We can also analyze the effect of moving in slightly 
from the edges of the phase diagram. Adding a small 
excess of ferromagnetic couplings in an XY spin-glass is 
not thought [16] to fundamentally affect the nature of 
the glass phase, so one does not expect to encounter any 
new phases moving in from the right edge of the phase 
diagram into the regime in which ^ > 0. Simi¬ 

larly, at non-zero T, the same is true near the left edge of 
the phase diagram where > 0 and T > 

We note that strictly speaking, in the framework of 
Eq. 5, the “globally s-wave” state exists only at non¬ 
vanishing temperatures. To see this, we consider the 
problem in terms of the Mattis transformed variables, 
so we consider the situation at low temperature, deep in 
the uniformly (ferromagnetically) ordered phase. Even 
though ^ there is some non-vanishing prob¬ 
ability (proportional to /J ) to find some pairs of 

nearby grains coupled by a strong negative (frustrating) 
Josephson coupling. Consequently, in the ground-state, 
there will be a small vortex loop enclosing this bond. It is 
easy to see that two such double-degenerate defects inter¬ 
act as dipoles, so they will in turn freeze into a dipolar 
glass phase at a temperature which is roughly propor¬ 
tional to their concentration, i.e. at least to first approx¬ 
imation Tgiass ^ Thus, at low enough tem¬ 
peratures and non-zero there occurs a phase with 

globally s-wave superconducting order coexisting with 
time-reversal symmetry breaking vortex glass order. At 
T > Tgiass the thermal fluctuations restore time reversal 
symmetry, leaving the globally s-wave state. 

Quantum fluctuations of the order parameter, which 
are not included in the model in Eq. 5, produce effects 
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FIG. 1. Schematic phase diagrams for a) = 0 

a nd J = + and b) = 0 and J = 

\/(+ (JG))2. The qualitative structure of the solid 
lines is justihed by asymptotic analysis in the text in the con¬ 
text of the (unproven) assumption that there is a hnite spin- 
glass transition temperature for the XY spin glass in d=3, 
while the dashed lines are included as a representative guess 
of how the phase diagram might be completed. The dotted 
lines represent the leading effect expected from quantum fluc¬ 
tuations - which are not explicitly included in the model. 


similar to non-zero T; quantum tunnelling between the 
two degenerate vortex loop states of an isolated defect 
thus effectively restore time reversal symmetry. This 
makes it possible for the globally s-wave state to exist for 
small enough even at T = 0 as indicated by the dot¬ 
ted line in Fig. la. It is not currently known how the su¬ 
perconducting phases merge at higher temperatures and 
intermediate and indeed it is likely non-universal. 

We have shown with dashed lines one possible, particu¬ 
larly simple completion involving a single tetra-critical 
point. 

In Fig. lb, we show a similar phase diagram for the 
situation in which = 0 while and are non¬ 
zero - this is relevant to the case in which the grains are 
all oriented, but have random shapes and separations. 
The nature of the phase diagram along the “edges” fol¬ 
lows from the same sort of analysis that led to Fig. la. 
Where the ^ we find the by now familiar glob¬ 
ally s-wave state, while for ^ we encounter the 
conventional d-wave state. At low temperatures, reflect¬ 
ing the intrinsic frustration when both and are 
comparable, both of these phases give way to a phase 
with coexisting glassy order (and hence subject to all the 
associated caveats). The phase diagram possesses a re¬ 
flection symmetry implied by the Mattis transformation 
in Eq. 7. Again, even the topology of the middle regions 
of the phase diagram is likely non-universal; with dashed 
lines we show a plausible minimal completion. The dot¬ 
ted lines indicate the region where (were we to include 
them in the model) quantum fluctuations of the order 
parameter destroy the glassy order of the vortex loops. 
Proximity effect and the Josephson couplings: 
We now address the origin of the various terms in 
Eq. 6. The Josephson coupling between grains of a d- 
wave superconductor separated by normal metal is con¬ 
trolled by the proximity effect. Andreev reflection at the 


superconductor-normal metal boundary generates a fi¬ 
nite value (0(^,7^)) inside the normal metal. We will 
consider the case where the grains are large compared to 
both the coherence length, and the mean free path 1. 

We begin by considering a single isolated grain em¬ 
bedded in a disordered normal metal matrix. It is im¬ 
portant to emphasize that, whether or not the order 
parameter inside the grain has d-wave character, the 
superconductor-metal interface breaks the point-group 
symmetry so a finite value of (0(r, r)) is induced in the 
neighboring metal. We will call the quantity (0(r, r)) 
the local s-wave component. It is this component which 
survives in the normal metal on distances larger than /, 
while the local d-wave and other-wave components decay 
exponentially at large distances. 

To obtain the requisite boundary condition for (0(r, r)) 
at the normal metal-superconductor interface, one has to 
match its values in the superconductor and in the normal 
metal close to the interface. In particular, its sign reflects 
the sign of (0(p±)) in the grain’s bulk for p = normal 
to the interface {\p\ = pi?). Thus, the sign of {(j){r^f)) 
changes along the boundary. 

The key consequence is illustrated graphically in Eig. 2, 
which is a schematic of a grain of a d-wave superconduc¬ 
tor embedded in a metallic matrix. At distances from 
the grain larger than I but small compared to the size of 
the grain, the sign of (0(r, r)) is dominated by the closest 
interface, and so it can be positive or negative depending 
on the orientation of that interface. However, far from 
the grain, the value of (0(r, r)) gets contributions from 
Cooper pairs that have scattered from all parts of the 
grain’s surface and then diffused to the observation point 
r. Thus its sign is determined by the sign of {(j){r^r))\s 
averaged over the entire surface of the grain. In other 
words at large distances, the sign of the anomalous aver¬ 
age is fully determined by the geometry of the grain. On 
the qualitative level the sign of (0(r, r)) is determined by 
the interference between the Cooper pairs traveling to 
the observation point r along different diffusive trajecto¬ 
ries, as illustrated by the two representative paths shown 
as the dashed lines in Eig. 2. 

To compute the Josephson coupling between a pair of 
grains formally requires the solution of this problem in 
the presence of two superconducting grains as a func¬ 
tion of their relative phases. However, it is straightfor¬ 
ward to understand the sign of the Josephson coupling by 
overlaying the patterns (^(r, r)) produced by each grain 
individually. As a result, for a pair of grains of compa¬ 
rable characteristic size R separated by a distance L, we 
find that for L > R, the sign of Rj is determined by a 
product single grain characteristics [12-14] - i.e. it ap¬ 
pears as a contribution to - while for L < R, its 
sign is determined by mutual aspects of the shape and 
orientation of the two grains - i.e. it appears as a con- 
tribution to Rp . Other than this, the considerations 
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controlling of the magnitude of the Josephson coupling 
are not substantially different from those in conventional 
S-N-S junctions. For instance, for L I 

Jij = ~ Vim ( 8 ) 

where D is the electron diffusion coefficient and = 
y^D/T the coherence length of the surrounding metal, 
and G is the characteristic conductance of the grain. In 
d = 3, generally G oc R so the term in parentheses scales 
as Jij ^ R? jL^. 

Importantly for present purposes, pairs of randomly 
oriented grains with R^ L contribute to while those 
with L ^ R contribute to Consequently, for fixed 

size grains, the material can be tuned across the phase 
diagram in Figs, la & b by varying the concentration of 
grains, as indicated. 

Identification of “globally s-wave” order: It 

should now be apparent that the global s-wave symme¬ 
try of the superconducting order refers to the symmetry 
of the induced order in the metallic host produced by 
the embedded d-wave grains. Even though the supercon¬ 
ductivity originates within the grains, the relative phase 
of order parameter from grain to grain is determined by 
the condition that the phase be constant throughout the 
metal. Another way to obtain an intuitive understanding 
is to imagine replacing the metal by an s-wave supercon¬ 
ductor. Now, consider the effect this has on the phase of 
the superconducting order in each (d-wave superconduct¬ 
ing) grain. It is clear that in this case, baring an acciden¬ 
tal degeneracy, the phase of the d-wave order parameter 
on each grain will be locked to that of the surrounding 
metal - either with the same phase or with a phase-shift 
of TT depending on the shape of the grain. If we now imag¬ 
ine continuously decreasing the strength of the intrinsic 
s-wave order in the matrix, by adiabatic continuity we 
would approach the situation we have discussed here. 
Further implications: In addition to its fundamental 
interest, the present results suggest new strategies for 
making better practical wires. Firstly, to avoid the vari¬ 
ous detrimental effects of frustration (and glassy phases), 
one would like to insure that this is ac¬ 

complished by insuring that the separation, L, between 
the grains is larger than or comparable to their char¬ 
acteristic radius, R. However, it is also desirable that 
the magnitude of the Josephson couplings be as large as 
possible. At the very least, this implies that we would 
like L < corresponding to the temperature (less than 
the bulk Tc) at which the wires are to be used. More¬ 
over, even when this inequality is satisfied, the coupling 
between two neighboring grains scales as J ^ R^ jL^. 
These various different considerations suggest that im¬ 
proved wires can be obtained by reducing the grain size 
and simultaneously increasing the concentration of grains 
subject to the condition L/R^l. In particular, there are 
likely regimes in which increasing L causes a transition 



Fig. 2 

FIG. 2. Schematic representation of the sign of the anoma¬ 
lous average of the pair creation operator produced in the sur¬ 
rounding disordered metal by the proximity effect coupling to 
a grain of a d-wave superconductor. The symbol inside the 
grain represents the structure of {(t>{p)) inside the grain. The 
dashed lines represent typical pair diffusion paths the con¬ 
tribute to the proximity effect. 

from the spin-glass to the globally s-wave regime, and 
correspondingly an increasing magnitude and decreasing 
time dependence of the critical current. 

Of course, this strategy has its limits - in order that 
the grains have undiminished local superconducting or¬ 
der, it is necessary that i? ^ and that R is large 
enough that quantum fluctuations of the superconduct¬ 
ing phase are negligible. Subject to this, for the cuprates, 
it would be particularly interesting to explore the situa¬ 
tion in which the London penetration depth. A, is large 
compared to the grain size, A > L ^ ^ in which 

case the superfluid density of the wire would be homo¬ 
geneous and isotropic, thus potentially mitigating some 
of the undesirable consequences of the quasi-2d nature of 
the cuprates. [20] 
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